EPJ manuscript No. 

(will be inserted by the editor) 



00 
^3 



Andreev Bound States in High Temperature Superconductors 

L. Alff, S. Kleefisch, U. Schoop, M. Zittartz, T. Kcmcn, T. Bauch, A. Marx, and R. Gross 
II. Physikalisches Institut, Universitat zu Koln, Ziilpicherstr. 77, D - 50937 Koln, Germany 



o 
o 

u 

Oh' 



o 



> 

o 
in 

o 

00 
On 
-I— > 



I 

O 

o 



Received: March 17, 1998 

Abstract. Andreev bound states at the surface of superconductors are expected for any pair potential show- 
ing a sign change in different fc-directions with their spectral weight depending on the relative orientation 
of the surface and the pair potential. We report on the observation of Andreev bound states in high tem- 
perature superconductors (HTS) employing tunneling spectroscopy on bicrystal grain boundary Josephson 
junctions (GBJs). The tunneling spectra were studied as a function of temperature and applied magnetic 
field. The tunneling spectra of GBJ formed by YBa2Cu307_5 (YBCO), Bi2Sr2CaCu208+i (BSCCO), and 
Lai.85Sro.i5Cu04 (LSCO) show a pronounced zero bias conductance peak that can be interpreted in terms 
of Andreev bound states at zero energy that are expected at the surface of HTS having a d-wave symmetry 
of the order parameter. In contrast, for the most likely s-wave HTS Ndi.85Ceo.i5Cu04-j; (NCCO) no zero 
bias conductance peak was observed. Applying a magnetic field results in a shift of spectral weight from 
zero to finite energy. This shift is found to depend nonlinearly on the applied magnetic field. Further con- 
sequences of the Andreev bound states are discussed and experimental evidence for anomalous Meissner 
currents is presented. 

Dedicated to J. Zittartz on the occasion of his 60th birthday 

PACS. 74.25.Fy Transport properties - 74.50.+r Proximity effects, weak links, tuimeling phenomena, and 
Josephson effects - 74.72.-h High-T sub c compounds - 74.76.Bz High-T sub c films 



1 Introduction 

Recent experimental and theoretical work on the symme- 
try of the order parameter in the high temperature su- 
perconductors (HTS) has led to the conclusion that for 
the majority of the cuprate superconductors the symme- 
try of the order parameter is dominated by a d^-i^y-i-wBMe 
component ||l|,^. This is in clear contrast to the metallic 
low-temperature superconductors for which the order pa- 
rameter is dominated by an isotropic s-wave component. 
The s-wave order parameter usually has the full symme- 
try of the point group of the underlying crystal structure 
and, hence, at the transition to the superconducting state 
only the global gauge symmetry is broken. Such pairing 
symmetry is denoted conventional. In contrast, the d-wave 
order parameter of the cuprate superconductors may not 
have the full symmetry of the point group of the crys- 
tal structure and, hence, beyond the global gauge sym- 
metry a further symmetry is broken at the transition to 
the superconducting state. Such pairing symmetry then 
is denoted unconventional. An important characteristic of 
the da;2_y2 -symmetry of the order parameter prevailing in 
most HTS is a sign change of the pair potential in orthog- 
onal /c-dircctions what is equivalent to a 7r-phase shift of 
the wave function. Usually the positive sign is taken along 
the a-axis of the unit cell and the negative sign along the 
6-axis. Furthermore, there are nodes of the pair potential 
in the [110] directions, where it changes sign. That is, the 



quasiparticles with corresponding fc-vector actually feel a 
vanishing pair potential what allows for low energy quasi- 
particle excitations. This is in contrast to a s-wave sym- 
metry, where a finite gap on the whole Fermi surface does 
not allow for quasiparticle excitations of arbitrarily small 
energy. We note that the actual order parameter present 
in hole doped HTS likely is formed by a mixture of a dom- 
inating (i3,2_j,2 -component and other components such as 
an s- or (ia;y-component. Furthermore, there is evidence 
that the order parameter of the the electron doped HTS 
Ndi.85Ceo.i5Cu04_j^ (NCCO) has a dominating s-wave 
component and no sign change in different fc-directions. 

There have been numerous experiments devoted to 
the determination of the symmetry of the order param- 
eter in the oxide superconductors. These experiments can 
be devided into those probing the amplitude by studying 
the quasiparticle excitation spectrum and those probing 
the phase of the order parameter in interferometer ex- 
periments employing multiply connected superconductors. 
Certainly, the key experiments have been the phase sensi- 
tive experiments designed by Tsuei and Kirtley which em- 
ploy HTS thin film tricrystals In these experiments 
a scanning SQUID-system has been used to show that 
there arc half-integer flux quanta in superconducting rings 
formed by three differently oriented HTS grains connected 
by grain boundary Josephson junctions. In this article we 
present a new class of phase sensitive experiments making 
use of the formation of Andreev bound states at surfaces 
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of HTS oriented parallel to the c-axis. An important con- 
sequence of an order parameter showing a sign change (or 
TT phase shift) in different fc-directions, is the formation of 
Andreev bound states at zero energy confined to the sur- 
face. Andreev bound states have been discussed first in the 
context of tunneling into unconventional superconductors 
by Buchholtz and Zwicknagel Q . Later on, the formation 
of Andreev bound states or midgap surface states having 
zero energy with respect to the Fermi energy and sizable 
areal density as a consequence of a (i3.2_y 2 -symmetry of 
the order parameter has been predicted by Hu These 
zero energy Andreev bound states can be probed by ab- 
plane tunneling spectroscopy and manifest themselves as 
a zero bias conductance peak. That is, the presence of 
Andreev bound states at zero energy represents definite 
evidence that the order parameter of the investigated su- 
perconductor changes sign along the Fermi surface. We 
emphasize that the same tt phase shift in the Josephson 
interference experiments by Tsuei et al. [Q is the origin of 
the Andreev bound states and the zero bias conductance 
peak in the a5-plane tunneling conductance. 



In this article we present comprehensive experimental 
data on a6-plane tunneling spectroscopy for various HTS 
materials. The data have been obtained using [001] tilt 
HTS grain boundary Josephson junctions (GBJs). In con- 
trast to low temperature scanning tunneling spectroscopy 
and experiments using SIN-type planar junctions, GBJs 
represent superconductor-insulator-superconductor (SIS) 
Josephson junctions |0,||j9j. An important advantage of 
the use of GBJs as compared to low temperature scan- 
ning tunneling spectroscopy is the very good long term 
stability of GBJs which allows for the detailed study of the 
temperature and magnetic field dependence of the tunnel- 
ing spectra. Furthermore, in GBJs one deals with inter- 
nal surfaces or interfaces. This implies that degradation 
or contamination effects due to ex situ processing in am- 
bient atmosphere are of minor importance. Only oxygen 
depiction in the region of the grain boundary can occur, 
but to a less extent as compared to bare surfaces used in 
low temperature scanning tunneling spectroscopy experi- 
ments. The tunneling data of the recent experiments using 
SIN-type junctions could be consistently explained by the 
formation of Andreev bound states as a consequence of 
a dominating d-wave symmetry of the order parameter 

~1,|T2|,B1. Here, we 



in the investigated HTS materials |10 



show that the same is true for the SIS-type GBJs fab- 
ricated from different HTS materials. Our experimental 
results can be compared to recent theoretical predictions 
by Barash et al. and Tanaka et al. on the Joseph- 
son behaviour of junctions with d-wave electrodes. Our 
detailed experimental study clearly shows that the zero 
bias conductance peak observed in HTS-GBJs is caused 
by Andreev bound states. Competing explanations for the 
origin of the zero bias conductance peak, in particular the 
magnetic scattering scenario, which is based on a model 
by Anderson and Appelbaum developed for NIN- 

type junctions containing magnetic impurity states, can 
be ruled out. 



Our analysis includes tunneling data obtained for the 
three hole doped HTS materials YBa2Cu307_a (YBCO), 
BiaSraCaCuaOg+s (BSCCO), and La2-a:Sr^Cu04 (LSCO). 
Here, YBCO has been investigated in the optimum doped 
phase with a critical temperature Tc of 90 K and in the 
underdoped phase with a Tc of about 60 K. Furthermore, 
the electron doped material NCCO has been studied. The 
analysis of this material is of particular interest with re- 
spect to the question whether or not there is a change 
in sign of the order parameter for NCCO. Up to now, 
there is convincing experimental evidence that NCCO has 
a dominating s-wave symmetry of the order parameter 
[0|l|,|2^,|l|, i.e., the order parameter of NCCO is ex- 
pected to show no change in sign for different fc-directions. 
Hence, Andreev bound states are expected only for hole 
doped HTS, which are believed to have a dominating d- 
wave symmetry of the order parameter, but should be ab- 
sent for NCCO. As we shall see below, this is indeed the 
case [^TB. 



2 Theoretical background 

There are several theoretical models that can account for 
a zero bias conductance peak in the tunneling spectra 
of SIN- or SIS-type junctions [||. Beyond the ABS sce- 
nario, the Appelbaum- Anderson model and the Blonder- 
Tinkham-Klapwijk model, phase diffusion in HTS Joseph- 
son junctions has been discussed |^,^. However, any 
conductance peak related to a supercurrent should be re- 
stricted to a much smaller voltage scale than the mV-scale 
usually measured for the zero bias conductance peak in 
most experiments. Moreover, any model based on super- 
currents fails for SIN-type junctions. In the following we 
therefore will restrict our discussion to the Appelbaum- 
Anderson model, , the Andreev bound state model, and 
the Blondcr-Tinkham-Klapwijk model. 



2.1 Appelbaum-Anderson model 

Zero-bias anomalies in the conductance of normal metal- 
insulator-normal metal (NIN) tunnel junctions often have 
been observed and can have various reasons |26,2^. In a 
model developed by Appelbaum and Anderson it is as- 
sumed, that the tunneling barrier is associated with local- 
ized paramagnetic states The tunneling conduc- 
tion electrons are exchange coupled to these states lead- 
ing to exchange scattering off the localized states. Similar 
to the s — d interaction model giving a resistance min- 
imum in dilute magnetic alloys (Kondo-type scattering 
[p8[), the Appelbaum-Anderson model can account for a 
zero bias conductance peak. The tunneling conductance 
G in the Appelbaum-Anderson model is given by three 
terms G = Gi + Ga + G3, where Gi is the contribution 
from all tunneling processes without spin interaction. G2 
is the spin exchange process contribution to the conduc- 
tance that, of course, depends on the applied magnetic 
field but not on the voltage at zero field. G3 is the Kondo- 
type contribution where an electron is scattered by the ex- 
change interaction leading to the interference of reflected 
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and transmitted waves. This additional tunneling channel 
contributes logarithmically to the conductance as 



G3(F,r) cxln 



En 



zV\ + nkBT 



(1) 



Here, Eq is an energy cut-ofF and n a factor close to unity. 
An applied magnetic field H results in a Zeeman splitting 
of the impurtiy states. Then, G2 is strongly suppressed 
for quasiparticle energies eV smaller than the distance 
between the Zeeman- levels and rises abruptly for \eV\ > 
gfj^sEl, where g is the Lande (/-factor, and the Bohr 
magneton. Furthermore, G3 is splitted by 26 — 2g^BH. 
Strictly speaking, G3 splits into three peaks with one peak 
centered at zero-bias. However, the spectral weight of this 
peak goes to zero as 6 is increased, probably impeding the 
measurement of this zero-energy peak . 

The Appelbaum- Anderson model has been successfully 
applied to Josephson tunnel junctions formed by low-tem- 
perature superconductors employing hydrogenated amor- 
phous silicon barriers Initially, it also was used to 
explain the experimentally observed zero bias conductance 
peak in tunnel junctions employing HTS electrodes Jl2| , ^0p i 
However, as will be shown below the Appelbaum- Anderson 
model cannot account for the temperature and magnetic 
field dependence of the zero bias conductance peak ob- 
served for the HTS junctions including the GBJs. We em- 
phasize that in general localized magnetic states in su- 
perconductors can lead to bound states within the gap as 
pointed out long ago by Zittartz et al. [|35| . 



2.2 Andreev bound state model 

The formation of Andreev bound states or midgap surface 
states in high temperature superconductors having zero 
energy with respect to the Fermi energy and sizable areal 
density as a consequence of a (i2;2_j^2-symmetry of the or- 
der parameter has been predicted by Hu The spectral 
weight of these bound states has a maximum for (110) 
oriented surfaces, whereas no such states are expected for 
(100) or (010) surfaces. That is, the mid gap states exist 
for any specular surface except for the lobe direction of 
the (ia;2_y2-gap perpendicular to the surface. The physical 
reason for the midgap states is the fact that quasiparti- 
cles reflecting from the surface experience a change in the 
sign of the order parameter along their classical trajectory 
and subsequently undergo Andreev reflection. Construc- 
tive interference of incident and Andreev reflected quasi- 
particles result in bound states confined to the surface. 
That is, the midgap states can be understood in terms of 
Andreev reflections, where the quasiparticle changes from 
par ticle-like to hole-like and vice versa, and k changes sign 
j^ , ^ . Therefore, the midgap states also are denoted as 
Andreev bound states. Reversal of the velocity by An- 
dreev reflection always is accompanied by a change in sign 
of the charge. Consequently, Andreev bound states can 
carry current, where current conservation is maintained 
by conversion of the bound state current to supercurrent 



far away from the surface. The Andreev bound states re- 
sult in a zero bias conductance peak in an a6-plane tunnel- 
ing spectrum. Thus, methods that probe the quasiparti- 
cle current on a surface or through an interface represent 
phase sensitive methods in the sense that the existence 
of Andreev bound states gives deflnite information on a 
change in sign of the pair potential. Hence, tunneling spec- 
troscopy becomes a valuable phase sensitive technique for 
probing the symmetry of the order parameter in HTS. 

The basic mechanism required to understand the for- 
mation of surface or interface bound states in dx-^-yi-WYS 
is the process of Andreev reflection . In general, An- 
dreev reflection always plays a role if there is a variation 
of the order parameter along the classical trajectory of a 
quasiparticle. That is, Andreev reflection can be viewed as 
scattering induced by spatial variations of the order pa- 
rameter. Such variations of the order parameter are well 
known for interfaces in artiflcial layer structures such as 
NS-interfaces. However, they also can occur along the tra- 
jectory of quasiparticles specularly reflected at the sur- 
face of a superconductor having a spatially anisotropic 
order parameter. In the process of Andreev reflection a 
|likc excitation undergoes branch conversion into 
ce excitation with reversed group velocity and vice 
versa. In order to discuss the importance of Andreev re- 
flection for the formation of Andreev bound states at sur- 
faces and interfaces of HTS in more detail let us consider 
the situation shown by Fig. |l|(a). An incident quasiparti- 
cle with a fc- vector corresponding to the negative lobe of 
the pair potential reflects specularly off the surface. For 
simpHcity, a (110) surface (a — 7r/4) and a quasiparticle 
that is incident under an angle /3 is considered. The order 
parameter Zi(k, r) experienced by the quasiparticle is a 
function of momentum k and space r. It is evident that 
incident and reflected wave packets propagate through dif- 
ferent order parameter fields. For the situation shown in 
Fig. |^(a) the order parameter has different sign for the 
incident and refiected wave packet. This is shown more 
clearly in Fig. |l](b), where we have plotted the order pa- 
rameter field along the classical trajectory of the quasi- 
particle. The rounding of A close to the surface is due to 
pair-breaking effects and plays no role for the formation of 
Andreev bound states. However, it becomes important if 
pair breaking frees significant spectral weight for subdom- 
inant pairing channels as will be discussed in section 2.3. 



In general, Andreev bound states occur at energies for 
which the phases of Andreev refiected particle-like and 
hole-like excitations interfere constructively. According to 
the Atiyah-Singer index theorem zero energy bound states 
are always formed, if the scattering induces a change in 
sign of the order parameter along the classical trajectory 
[ p8| . For the situation shown in Fig. |l| this is always di- 
rectly at the interface (d = 0). Thus, changing the incident 
angle (3 will not change the location of the Andreev bound 
states at d = 0. Reducing a from 7r/4 towards a = re- 
sults in a increasingly smaller range of incident angles (3 
for which a change in sign of the order parameter field 
is obtained. Thus, the spectral weight of the zero energy 
bound states decreases with decreasing a and vanishes for 
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Fig. 1. Formation of Andreev bound states: In (a) a mixed 
real/momentum space representation is shown for a (110) ori- 
ented surface of a HTS having a d^2_y2 symmetry of the order 
parameter. The propagation direction of the quasiparticle is 
indicated by the arrows and can be transferred to a direction 
in fc-space by plotting a k- vector from the origin of the lobe to 
its edge parallel to the propagation direction in real space. In 
(b) the order parameter A is plotted versus the path d along 
the quasiparticle trajectory. For the case of a sign change on 
reflection the Andreev bound states are always formed at the 
surface {d — 0). In (c) a schematic illustration of the differ- 
ent reflection and transmission processes of quasiparticles at a 
grain boundary (GB) interface is given. 



a = 0, i.e. for a (100) or (010) oriented surface. That is, 
the spectral weight of the Andreev bound states is largest 
for (110) surfaces and decreases continously towards zero 
for (100) or (010) surfaces. It is evident that for trajec- 
tories orthogonal to the Cu02-plane no Andreev bound 
states are formed. 



Above we have considered the situation present at the 
surface of a HTS. We now briefly discuss the situation for 
a SIS-structure present for the GBJs studied in our ex- 
periments. In Fig. |l|(c) the possible transmission and re- 
flection processes at a grain boundary interface are shown 
schematically. An incident particle-like excitation is in- 
jected from the left bulk d-wave HTS under a flnite angle 
with respect to the grain boundary plane. Discussing the 
reflection processes, the particle-like excitation can be ei- 
ther specularly or Andreev reflected. In the latter case it 
is turned into a hole-like excitation propagating in oppo- 
site direction. In the same way, transmission through the 
insulating grain boundary barrier can yield particle- and 
hole-like excitations. As described above for the case of 
reflection at a surface, the quasiparticles involved in the 
different reflection and transmission processes experience 
different pair potentials depending on their fc-vector and 
the relative orientation of the order parameter in both 
electrodes with respect to the grain boundary barrier in 
both electrodes. We also note that Cooper-pairs have to 
be involved to establish charge conservation. 

We point out that different order parameter symme- 
tries will produce different fc-dependencies of the spec- 
tral weight of Andreev bound states. Of course, Andreev 
bound states are present only if there is a change in sign 
of the pair potential in different fc-directions and are ab- 
sent if there is none. For example, for a d^y- and d^2_y2- 
symmetry the maximum spectral weight of the Andreev 
bound states is obtained for a (100) or (110) oriented sur- 
face or interface, respectively. In contrast, for a s-wave 
symmetry no sign change of the order parameter and, 
hence, no Andreev bound states are present. Thus, the 
experimental observation of Andreev bound states gives 
deflnite evidence for an order parameter changing sign on 
the Fermi surface and represents a phase sensitive probing 
technique. In cases where the order parameter symmetry is 
not fully established, as for example for the electron doped 
NCCO, Andreev bound states sensitive experiments can 
be used to clarify whether or not the order parameter 
changes sign in different fc-directions. 



2.3 Consequences of the Andreev bound state model 

There are several new and interesting phenomena that are 
directly related to the existence of Andreev bound states 
at surfaces and interfaces of HTS that will be addressed 
in the following. So far, there is significant experimental 
evidence only for the zero bias conductance peak in tun- 
neling data due to Andreev bound states. However, other 
effects like an anomalous temperature dependence of the 
critical current of HTS Josephson junctions or an anoma- 
lous temperature dependence of the London penetration 
depth due to "anti-Meissner" currents have not yet been 
experimentally confirmed. 



2.3.1 Zero bias conductance peak 

Since Andreev bound states can carry current the most 
evident experimental consequence is that they should pro- 
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duce a pronounced zero bias conductance peak in the ab- 
plane tunneling conduction of experiments involving at 
least one d-wave HTS electrode. For the interpretation of 
the conductance of superconductor-insulator-normal metal 
(SIN) junctions that can be formed by using a normal 
metal tip in low temperature scanning tunneling spec- 
troscopy, several groups extended the calculations by Hu 
1^. In particular, expressions for the density of states and 
the conductance of such junctions employing d-wave HTS 



have been derived [|1^ 
nent feature in the ca. 



15[|39j,|^,|l]j4|]. The most promi- 
culated conductance versus volt- 
age curves indeed was a pronounced zero bias conduc- 
tance peak as a direct consequence of the Andreev bound 
states. Recently, these states also have been observed ex- 
perimentally using low temperature scanning tunneling 
spectroscopy JT0p3p^ , in particular on well defined (110) 
oriented HTS surfaces |11 . We note that Andreev bound 
states are absent at surfaces perpendicular to the c-axis di- 
rection. Therefore, no zero bias conductance peak has been 
observed by Renner et al. in low temperature scanning 
tunneling spectrosco py exp eriments performed on (001) 
oriented surfaces |4^^,^. A further possibility for the 
realization of SIN-junctions is the formation of planar 
junctions using a normal metal and a HTS electrode sep- 
arated by a thin insulating layer. Also for such junctions a 
zero bias condvictance peak has been observed in most ex- 

i|,|^,|i|]3: 



periments |12 



Up to now the theoretically 



predicted zero bias conductance peak has been found for 
both SIN- and SIS-type junctions ||,|9|j5|,|l| . Below we 
present comprehensive experimental data on the tunneling 
conductance of various high temperature superconducting 
GB Js giving a very complete view of the zero bias conduc- 
tance peak as a consequence of Andreev bound states. We 
note that non zero-energy states are also possible as a con- 
sequence of a suppression of the order parameter at the 



surface |42 



2.3.2 Surface roughness 

In our discussion of the Andreev bound state model we 
always have assumed perfectly flat surfaces and interfaces 
that usually are not present in any experimental situa- 
tion. The spectral weight of the Andreev bound states and, 
hence, the measured zero bias conductance peak is very 
sensitive to surface roughness on a length scale larger than 
the coherence length. Since the coherence length of HTS 
is of the order of 1 nm, almost all surfaces can be consid- 
ered rough. The influence of surface roughness has been 
theoretically calculated by several authors [^ , ^2l|5^ , |5^ . 
The basic consequence of surface roughness is that one 
has to average over different surface orientations. As a 
consequence, for a (i3;2_y 2 -symmetry of the order param- 
eter and a (110) oriented surface, for which the spectral 
weight of the Andreev bound states is maximum, the zero 
bias conductance peak is suppressed by surface roughness. 
Vice versa, a finite zero bias conductance peak can appear 
for (100), (010), and (001) surfaces and interfaces due to 
surface roughness. That is, surface roughness smears out 
the pronounced dependence of the Andreev bound states 












Fig. 2. Qualitative plot of the amplitude of the order parame- 
ter as a function of the distance from the surface for a (100) (a) 
and (110) (b) surface without (i) and with (ii) surface rough- 
ness, respectively. Also shown is a surface induced s-wave order 
parameter that can be obtained for a (110) surface by surface 
pair breaking effects for the dominant d-wave and subsequent 
pairing into the subdominarit s-wave channel. 



spectral weight on the relative orientation of the surface 
and the order parameter. As an example, in Fig. ^ the 
order parameter is plotted versus the distance from the 
surface for smooth and rough surfaces with different ori- 
entation. It is evident that due to surface roughness differ- 
ently oriented surfaces behave in a similar way. This makes 
it experimentally more difficult to distinguish the differ- 
ent directions within the superconducting plane. However, 
since the sign change of the pair potential is a necessary 
prerequisite for the formation of Andreev bound states, 
any amount of surface roughness can not produce a zero 
bias conductance peak for an order parameter having no 
sign change on the Fermi surface such as an s-wave order 
parameter. 



2.3.3 Anomalous temperature dependence of the Josephson 
current density 

A further important consequence of the Andreev bound 
states is an anomalous behaviour of the temperature de- 
pendence of the maximum Josephson current Ic at low 
temperatures ||55|,|56| . Recently, is was predicted by Tanaka 
and Kashiwaya |55|| that for a fixed phase difference be- 
tween the electrodes of a Josephson junction formed by 
d-wave superconductors the supercurrent across the junc- 
tion can be either positive or negative depending on the 
injection angle of the quasiparicles. Loosly speaking, the 
junction can be thought as a combination of so-called 
0- and TT-junctions. Due to different temperature depen- 
dencies of both components of the critical current, this 
may result in a nonmonotonous temperature dependence 
of the Josephson current. In particular, such behaviour 
is expected for symmetric [001] tilt HTS grain boundary 
junctions. We point out that this feature is expected only 
when Andreev bound states are formed at the interface of 
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both junction electrodes and if there is a change in sign of 
the Josephson current as a function of injection angle. The 
predicted nonmonotonous Ic{T) is in clear contrast to the 
behaviour of traditional Josephson junctions formed by s- 
wave superconductor, for which a monotonous increase of 
Ic with decreasing temperature is theoretically predicted 
and experimentally observed. The nonmonotonous Ic(T) 
dependence for d-wave Josephson junctions again is sen- 
sitive to interface roughness. For increasing roughness the 
anomaly is shifted to lower temperatures [ p6| . This may be 
the reason why this effect has not yet been convincingly 
experimentally observed. Additionally, the presence of a 
subdominant s-wave component of the order parameter 
as discussed in the next section is predicted to suppress 
strongly the enhancement of the Josephson current at low 
temperatures [p7[. 



2.3.4 Splitting of the zero bias conductance peak by a 
broken time reversal symmetry state 

It has been shown in general that a locally time reversal 
symmetry breaking state at interfaces of unconventional 
superconductors can exist that decays exponentially to- 
ward the bulk In a broken time reversal symmetry 
state the Andreev bound states shift to finite energy re- 
sulting in a split zero bias conductance peak. A broken 
time reversal symmetry state is obtained by applying a 
magnetic field. However, even at zero applied magnetic 
field time reversal symmetry can be broken by the pres- 
ence of two order parameters at the surface having a 7r/2 
relative phase difference. The possible coexistence of a sec- 
ond order parameter with s-wave symmetry that can be- 
come stable at surfaces of HTS but is completely dom- 
inated by the d-wave potential in the bulk material has 
been discussed by Sigrist et al. |58 and Matsumoto et al. 
P^ . Such situation is possible at the surface of a d-wave 
HTS, since the surface may act as strong pair breaker for 
the d-wave paired quasiparticles thereby releasing spectral 
weight for the formation of pairs in a subdominant (e.g. 
s-wave) pairing channel. Below a surface transition tem- 
perature then a surface order parameter can develop that 
spontaneously breaks time reversal symmetry and results 
in spontaneous surface currents. The pairing symmetry at 
the surface can be a mixed d + is (or in another nota- 
tion Big + iA2g) phase. Note that both order parameters 
coexist at the surface with a relative phase shift of 7r/2 
that leads to a spontaneous surface current. The Andreev 
bound states then are shifted to finite energy due to the 
Doppler shift • Ps , where Ps is the condensate momen- 
tum due to the spontaneous supercurrent. This in turn 
leads to a split zero bias conductance peak in zero applied 
magnetic field [ ^9|j6C|] . That is, the measurement of the 
conductance in a5-plane tunneling experiments can give 
information on the presence of a surface order parameter 
in HTS that spontaneously breaks time reversal symme- 
try. Recently, this effect has been found to be consistent 
with the t — J-model in the case of suffi ciently large 
J at the surface of a d-wave superconductor 



Recently, a zero field splitting and a nonlinear evo- 
lution of the magnitude of the splitting with increasing 
applied magnetic field as predicted by Fogelstrom et al. 

has been reported for planar SIN-junctions [p^^ 
and for low temperature scanning tunneling spectroscopy 
experiments pJ. 



2.3.5 Anomalous Meissner current 

The London-penetration depth Al of a superconductor is 
known to be related to the density of paired quasiparticles. 
For the different symmetries of the order parameter a dif- 
ferent behaviour of A^fT) has been predicted and experi- 
mentally found pq , |66|j67|| . The presence of Andreev bound 
states at surfaces of HTS can have a significant effect on 
the low temperature dependence of the London penetra- 
tion depth. In order to clarify this effect let us consider the 
surface of a d-wave superconductor where Andreev bound 
states are formed. In a state with established time reversal 
symmetry there is no net current parallel to the surface, 
since in thermal equilibrium the contributions of all An- 
dreev bound states exactly cancel each other. However, 
applying a magnetic field parallel to the surface results 
in a broken time reversal symmetry state and a Meissner 
shielding current. This in turn results in a Doppler shift 
V/ • , where Ps now is the condensate momentum due to 
the Meissner shielding current. Since the Doppler shift has 
opposite sign for Andreev bound states with opposite Vf, 
in thermal equilibrium there is no longer an equal popu- 
lation of states with opposite Vf. This results in a net sur- 
face current that is opposite to the usual Meissner shield- 
ing current and therefore is denoted as "anti-Meissner" 
current. In contrast to the usual Meissner shielding cur- 
rent flowing in a surface layer of depth Al, this "anti- 
Meissner" current flows in a surface layer given by the ah- 
plane coherence length £^ab- It becomes important at low 
temperatures where is causes an anomalous temperature 
dependence of the London penetration depth \l{T) [Q. 
Due to the "anti-Meissner" current the London penetra- 
tion depth is slightly increased beyond the value expected 
without this current. First experimental evidence for this 
effect has been reported recently by Walter et al. |69| and 
more evidence will be presented below for HTS-GBJs. 



2.4 Blonder-Tinkham-Klapwijk model 

Blonder, Tinkham, and Klapwijk (BTK) developed a gen- 
eralized Andreev reflection model for junctions involving 
superconductor-normal metal (S/N) interfaces [^0|j7l|. To 
account for a finite barrier strength at the S/N interface 
they introduced a dimensionless factor Z . For example, in 
the limiting case of Z = the transmission coefficient of 
the junction that is given as + Z'^) becomes unity. 
In the Z — limit and for small values of Z the BTK 
model predicts a peak in the normalized differential con- 
ductance at zero voltage that reaches a maximum value of 
2. This peak decreases monotonously with increasing volt- 
age, that is, the BTK model predicts a conductance peak 
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at zero bias. However, no gap structure in the conduc- 
tance versus voltage curves is expected to be observed for 
small Z. That is, within the BTK-model the observation of 
both a gap like structure and a zero bias conductance peak 
is not expected. We emphasize that the BTK-model was 
developed for conventional BCS superconductors. There- 
fore, in its original form the BTK-model is not sufficient 
to explain the experimentally observed tunneling spectra 
of d-wave superconductors. However, the BTK-model can 
be used as a starting point for more elaborate theories. 
The Andreev bound state model represents a generalized 
BTK-type model that accounts for an arbitrary Z and an 
anisotropic pair potential. 



3 Experimental Results 

3.1 Sample preparation and experimental techniques 

In our experiments both symmetric (ai = — 012) and asym- 
metric (ai =0; 02) [001] tilt HTS-GBJs fabricated on 
SrTiOs bicrystal substrates were used. A sketch of the 
grain boundary junction geometry is shown in Fig. |[ The 
symmetric GBJs had a total misorientation angle (ai -1-0:2) 
of 24° and 36.8°, for the asymmetric [001] tilt GBJs we 
used a-2, = 45°. Only c-axis oriented epitaxial thin films 
were grown on the bicrystal substrates. In this way the 
tunneling direction in the grain boundary junction con- 
figurations always is along the afe-plane. The YBCO thin 
films were deposited using both hollow cathode magnetron 
sputtering and pulsed laser deposition from stoichiomet- 
ric targets. YBa2Cu307_5 samples were fabricated with 
different oxygen content. Both fully oxidized YBCO with 
a critical temperature Tc of about 90 K and oxygen defi- 
cient YBCO with Tc values ranging between 55 and 60 K 
were used. We note that the reduction in the oxygen 
deficient YBCO is obtained by well controlled oxygen de- 
pletion and is not caused by the presence of magnetic 
impurities in the samples or at the grain boundary. Fur- 
thermore, BSCCO-GBJs were fabricated in the same way 
with a Tc of the BSCCO-films of about 80K |7|. This 
Tc value is reduced by about 10 K as compared to opti- 
mum doped samples The epitaxial LSCO thin films 
have been grown on the bicrystal substrates by reactive 
coevaporation from metal sources using ozone as reaction 
gas |7^,^,Q . The Tc of these samples typically was about 
24 K. This value is comparable to other epitaxial thin film 
samples but slightly reduced compared to that of the best 
single crystals of LSCO. This most likely is due to strain at 
the interface between the substrate and the epitaxial thin 
film [|7^,|7^. Finally, the NCCO films were deposited simi- 
lar to the LSCO films by molecular beam epitaxy (MBE) 
using ozone as oxidation gas [Q. The Tc of about 24 K 
is one of the highest reported for this material in litera- 
ture including even NCCO single crystals [Q. After their 
growth, the HTS thin film bicrystals have been patterned 
using optical lithography and Ar ion beam etching to ob- 
tain a GBJ structure as sketched in Fig. |. 

The measurements of the current-voltage [/(V^)] and 
conductance vs. voltage [0(1^) = dl{y) j dV\ characteris- 




Fig. 3. Sketch of the sample geometry of a symmetric [001] 
tilt HTS-GBJ. The crystal axis a, b, and c and the misorien- 
tation angles qi and a2 are indicated. For a symmetric and 
asymmetric junction we have a\ = —02 and |ai| 7^ |a2], re- 
spectively. Also given is a magnified view of the grain bound- 
ary microstructure indicating the strong faceting of the grain 
boundary. The nominal misorientation angle is only the aver- 
age angle of a wide range of different angles of the individual 
facets. 



tics were performed in a standard four-probe arrangement. 
If necessary the critical current was suppressed by apply- 
ing a small magnetic field parallel to the grain boundary 
plane corresponding to a minimum in 1^ vs. H . Zero field 
measurements were performed in magnetically shielded 
cryostats placed in an rf-shielded room. Also, high mag- 
netic fields up to 12 T and low temperatures down to 
100 mK were used in our study. 



3.2 Transport mechanism in GBJs 

The electrical transport properties of grain boundaries in 
the cuprate superconductors have been studied intensively 
over the last years. Recent reviews have been given by 
Gross et al. jQ,^. Here, we briefly summarize the main re- 
sults relevant to the tunneling spectroscopy experiments 
discussed below. Phenomenologically, the superconduct- 
ing properties of the GBJs can be well modeled by the re- 
sistively and capacitively shunted junction (RCSJ) model 
[0IZll ■ In order to describe the vaste majority of the present 
experimental data on the electrical transport and noise 
properites of HTS-GBJs the intrinsically shunted junction 
(IS J) model has been suggested l ]80| , ^ . In this model it 
is assumed that there is a thin insulating grain boundary 
barrier containing a large density of localized defect states. 
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The microscopic transport mechanism then is direct and 
resonant tunneling via the locahzed states. A key feature 
of the ISJ-modcl is the fact that the transport of Cooper 
pairs via the locahzed states is prevented by Coulomb re- 
pulsion, i.e. Cooper pairs have to use the direct tunneling 
channel, whereas due to the large density of defect states 
the quasiparticle transport is dominated by resonant tun- 
neling via a single localized state [[73|j8^j8^ , p^ , ^5l|8^ , |87| . 
The resonant channel can be viewed as to provide an in- 
trinsic resistive shunt leading to the name intrinsically 
shunted junction. The ISJ-model well explains the scaling 
behaviour Vc = JcPn {JcY with q ^ 0.5 observed for 
the GBJs ||80| , |88| . Here, Jc and p„ are the critical current 
density and the normal resistance times area of the GBJ. 
Furthermore, it naturally accounts for the large amount 
of low frequency excess noise due to trapping and release 
of charge carriers within the localized states | |7^ , |86| . 

We emphasize that the tunneling processes present in 
GBJs are elastic processes. If inelastic transport processes 
would be important, a significant temperature dependence 
of the tunneling conductance for voltages well above the 
gap voltage would be expected. However, the experimental 
data on GBJs show a temperature independent conduc- 
tance |]3^ ; i. e. inelastic tunneling via two and more local- 
ized centers is negligible. Hence, the tunneling quasipar- 
ticles can carry spectroscopic information on the density 
of states (DOS) in the junction electrodes and thereby 
e.g. on the energy gap. In contrast, such spectroscopic 
information is destroyed, if multistep inelastic tunneling 
or even variable range hopping would be the dominating 
transport process. We also note that the localized states 
within the grain boundary barrier have an about white 
energy distribution. Hence, they do not act as an energy 
filter producing structures in the conductance versus volt- 
age curves of HTS-GBJs. Because all junctions used in our 
experiments could be well fitted within the above picture, 
we assume that the same tunneling mechanism is present 
in all investigated junctions. 

There have been extensive studies on the microstruc- 
ture of the grain boundary barrier. High resolution trans- 
mission electron microscopy (HRTEM) revealed a strongly 
faceted grain boundary barrier with facets on a nanoscale 
p9| . This means that the macroscopic misorientation an- 
gles Qfi and a2 of the grain boundary, which are deter- 
mined by the bicrystal substrate, represent average values 
only. However, on a microscopic scale these angles can vary 
considerably due to the faceting with \ai \ + \a2\ = const. 
As a consequence in some cases GBJs formed by c?-wave 
HTS materials have to be modelled as a parallel array 
of so-called "0"- and "tt" -junctions showing an irregular 
magnetic field dependence of the critical current | |9C| , |9H . 
This effect is negligible for symmetric GBJs with small 
misorientation angles (|q;i| -I- \a2\ < 25°) however, 
it becomes more pronounced for increasing misorienta- 
tion angle and stronger faceting. In general, the faceted 
grain boundary barrier can be viewed as a rough interface. 
The impact of surface roughness on the spectral weight of 
Andreev bound states on their consequences already has 
been discussed above. Unfortunately, up to now no exact 




V (mV) 

Fig. 4. (a) Normalized differential conductance vs. voltage 
for a 24° [001] tilt YBCO-GBJ for temperatures between 4 
and 51 K. The sample was underdoped with Tc « 55 K. Some 
data points close to zero are removed because of the Joseph- 
son current peak, (b) The same dependence for temperatures 
between 45 and 51 K on an enlarged scale. One clearly can 
distinguish the sharp Josephson current peak at zero bias, the 
much broader zero bias conductance peak, and the supercon- 
ducting gap structure even for T close to T^. 

quantitative description of the effects of a faceted grain 
boundary on the electrical transport properties is avail- 
able. Below we will qualitatively discuss possible effects 
of the faceting that can be understood from simple angle 
averaging. 



3.3 Tunneling spectroscopy of HTS-GBJs 

In this section we present our experimental results on the 
tunneling spectroscopy on various HTS materials using 
GBJs. In the next section we will discuss in detail our 
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results in the framework of the different models discussed 
above. 

We first consider the conductance versus voltage curves 
of the hole doped HTS YBCO, BSCCO, and LSCO, which 
are believed to have a d-wave symmetry of the order pa- 
rameter. Fig. IJa shows a set of typical G'(y)-curves nor- 
malized to the normal state conductance G„. The curves 
were obtained for a 24° [001] tilt YBCO-GBJ at tempera- 
tures ranging between 4 and 51 K. The Joscphson current 
peak at zero voltage has been removed for clarity. The 
sample of Fig. |a was underdoped and had a Tc of about 
55 K. In Fig. ^b some spectra are shown on an enlarged 
scale in order to better show the presence of the broad zero 
bias conductance peak up to the Tc of the material. Fur- 
thermore, the superconducting gap structure and a nar- 
row peak at y = due to the Josephson current can be 
observed up to Tc- For other hole doped HTS materials 
such as BSCCO and LSCO the same behaviour was ob- 
served. For LSCO with Tc ~ 24K the experimental results 
have been published recently . For all samples the zero 
bias conductance peak is found to decrease with increas- 
ing temperature and to disappear at Tc. Furthermore, also 
the gap structure is found to disappear at Tc- 

The measured tunneling spectra are almost perfectly 
symmetric about zero voltage as expected for a symmet- 
rical tunnel junction with electrodes formed by the same 
material. Far above the gap voltage the conductance was 
found to be independent on temperature and has an about 
parabolic shape pq] . The parabolic shape can be explained 
by the effect of the applied voltage on the shape of the bar- 
rier potential. The background conductance can be viewed 
as normal state conductance Gn ■ In order to show the tem- 
perature dependence of the spectra more clearly, the con- 
ductance data have been normalized with respect to G„. 
Fig. ^ clearly shows that below the gap voltage the den- 
sity of states is reduced by about 30%. With increasing 
temperature the conductance smoothly merges towards 
to the normal state conductance when the temperature 
is approaching Tc. Recently, Ekin et al. |^ reported the 
disappearance of the zero bias conductance peak at tem- 
peratures well below the transition temperature of the ma- 
terial. Such behaviour has not been observed in our study 
and thus might be an effect of sample fabrication. In our 
case a clear correlation between the appearance of a mea- 
surable zero bias conductance peak and Tc was observed. 

The materials we used for our experiments were both 
hole doped HTS (60 and 90 K phase of YBCO, BSCCO, 
and LSCO) and an electron doped HTS (NCCO). We now 
compare the data of the hole doped HTS to that obtained 
for the electron doped HTS. There is significant exper- 
imental evidence that the hole doped materials have a 
c?-wave symmetry of the order parameter, whereas NCCO 
has a dominating s-wave component. From our experi- 
ments on GBJs fabricated from these materials we can 
conclude that all hole doped HTS show a pronounced zero 
bias conductance peak, whereas such peak is completely 
absent for the electron doped NCCO. In order to demon- 
strate this observation, in Fig. ^ we show data for fully 
oxidized YBCO, as well as for LSCO and NCCO. Again 




Fig. 5. Normalized conductance vs. normalized voltage of 
GBJs formed by LSCO, YBCO (90 K phase), and NCCO. 
The YBCO-GBJ was an asymmetric [001] tilt GBJ while the 
LSCO- and NCCO-GBJ were symmetric [001] tilt GBJs. For 
the NCCO-GBJ the Josephson current peak has been removed. 



the data is normalized to the normal conductance G„. 
Furthermore, the voltage scale is normalized to the gap 
voltage Vg in order to compare data of materials with dif- 
ferent gap values. Here, we have chosen Vg = Aq/c and 
not, as for tunnel junction using conventional BCS-like su- 
perconductors, Vg = 2Z\o/e. We also note that the exact 
position of the conductance peak due to the gap depends 
on the orientation of the electrodes with respect to the 
barrier. This is a consequence of the (interface) density 
of states in a d-wave superconductor |14] . In general, the 
same consideration should apply for a highly anisotropic 
s-wave superconductor. Therefore, in a first order approx- 
imation Vg ~ zio/e is assumed for NCCO as well. Typical 
gap voltages of the different materials were 6 mV for both 
LSCO and NCCO, 15 mV for the 60 K phase of YBCO 
(YBCO-60), 20 mV for the 90 K phase (YBCO-90), and 
25 mV for BSCCO. 

With respect to the gap structure in the measured tun- 
neling spectra we would like to add an interesting obser- 
vation concerning the magnitude of the gap. As can be 
seen from Fig. ^ the magnitude of the gap seems to be 
about constant with increasing temperature. However, the 
states within the gap are filled up with increasing tem- 
perature until the normal state value is reached. Above 
Tc the presence of a gap could not be observed in our 
experiments. A similar behaviour has been observed by 
Renner al. |^ , ^ using low temperature scanning tunnel- 
ing spectroscopy. However, in contrast to our observation 
these authors report a gap persisting even above Tc. We 
note, that the tunneling direction in the experiments by 
Renner et al. was along the c-axis, whereas it is along 
the a6-plane in our study. The observed temperature evo- 
lution of the gap structure completely deviates from the 
BCS-behaviour and is probably related to the pseudogap 
that has been observed in the HTS ||3|j94). 

We next consider the magnitude of the gap as a func- 
tion of the doping level. Recently, an increasing gap was 
found with decreasing doping concentration for BSCCO 
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1^,^. In the case of YBa2Cu307_5 we investigated the 
two doping concentrations (5 ~ and S ~ 0.4. In contrast 
to the recent work on BSCCO, we observed a hnear scal- 
ing of the gap values with Tc- We note, however, that this 
behaviour has to be examined in more detail in future. 
Also for NCCO, the gap value was found to scale with Tc- 
A junction with a Tc of about 14 K had a gap value of 
about 3.5nieV, whereas junctions with ~ 24K showed 
gap values of about 6meV. 

The conductance versus voltage curves of GBJs fabri- 
cated from the c?-wave HTS |Q all show some common 
properties. Firstly, they have a pronounced zero bias con- 
ductance peak. Secondly, the width of the zero bias con- 
ductance peak is always about several mV and is not nar- 
rowing with decreasing temperature. Thirdly, the conduc- 
tance is reduced in an energy range corresponding to the 
superconducting gap energy most likely due to a reduced 
density of states. For all d-wave HTS only a reduction of 
the density of states of about 30 to 50% is observed. Ad- 
ditionally, a broad and flat peak in the G{V) curves is ob- 
served at the gap voltage, however, not for all samples. Fi- 
nally, reducing the temperature to very small values (down 
to 100 mK) results in a zero bias conductance peak with a 
shape that is more Lorentzian-like than Gaussian-like. As 
shown in Fig. |^, a Gaussian shape is obtained for LSCO at 
4.2 K, whereas for YBCO-90 already a Lorentzian shape 
is obtained at this temperature. This indicates that the 
measured linewidth is not determined by extrinsic effects 
such as thermal smearing or the finite experimental res- 
olution but represents an intrinsic linewidth of the zero 
bias conductance peak. Galculating the linewidth of the 
zero bias conductance peak, an imaginary part can be 
added to the energy introducing a phcnomcnological life- 
time parameter F of the order of Ao/10 |l^. Recently, a 
momentum-dependent intrinsic broadening of the surface 
bound states has been predicted |p7|| . 

The behaviour of NGGO-GBJs differs from that of 
GBJs fabricated from the d-wave HTS in several features. 
Firstly, the G{V) curves of optimum doped NCCO-GBJs 
with a maximum Tc of 24 K show a much clearer peak 
structure at Vg than those of the d-wave HTS-GBJs. Sec- 
ondly, the gap structure is more pronounced and the den- 
sity of states is reduced to about 40% of the normal state 
value. For a NCCO-GBJ with a reduced Tc both the peak 
at Vg and the reduction of the density of states becomes 
smaller. However, up to now no detailed tunneling spec- 
troscopy as a function of doping has been performed. Fi- 
nally, it is interesting to note that the spectra of NGGO 
have, apart from the zero bias conductance peak, more 
similarities to those of the hole doped HTS materials than 
to those of conventional BCS-superconductors. 

We finally consider the magnetic field dependence of 
the conductance versus voltage curves. The magnetic field 
dependence of the tunneling spectra of a YBGO-GBJ (60 
K phase) is shown in Fig. ^(a). The magnetic field H was 
applied parallel to the grain boundary barrier in c-axis di- 
rection. The magnetic field was increased in 0.5 T steps. In 
Fig. ^(a) the Josephson current has not been suppressed 
explaining the sharp spike at zero voltage in the zero field 
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Fig. 6. (a) Magnetic field dependence of the zero bias con- 
ductance peak at 4.2 K for a symmetric [001] tilt YBCO-GBJ 
(60 K phase). In (b) {G{V, H) - G(l/, 0) is plotted for the same 
sample. The applied magnetic fields ranged between and 12 T 
(0.5 T steps). 

G{V) curve on top of the much broader zero bias conduc- 
tance peak caused by the Andreev bound states. The effect 
of the applied magnetic field is the suppression of spectral 
weight around zero energy. The spectral weight removed 
around zero energy is shifted to higher energies about 
the whole energy range. Hence, in the G{V, H) — G{V, 0) 
vs. voltage curves shown in Fig. |^(b) two peaks appear 
with the peak position increasing non-linearly with in- 
creasing applied magnetic field |^ . We note that the over- 
all spectral weight at 12 T is reduced by about 3% below 
the value at zero magnetic field if we consider a voltage 
interval of ±12 mV. However, it is not likely that the con- 
servation of the number of states is violated but that the 
missing spectral weight is shifted to even higher energies 
not contained in the considered voltage interval. For LSCO 
at 100 mK qualitatively the same result has been obtained 
[ p2[ . Note that the thermal smearing is only of the order 
of 10 fieV at 100 mK. This explains that for LSCO the 
zero bias conductance peak in the conductance vs. volt- 
age curves changes from a Gaussian shape at 4.2 K to 
Lorentzian shape at 100 mK. In contrast to reference 
but in accordance to reference |^ no direct splitting of 
the zero bias conductance peak was found at any applied 
magnetic field up to 12 T. 
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3.4 Comparison to low temperature scanning 
tunneling spectroscopy measurements 



In this subsection we briefly compare the results obtained 
for SIS-type HTS-GBJs to those obtained using low tem- 
perature scanning tunneling spectroscopy, where a SIN- 
configuration is used. A low temperature scanning tun- 
neling spectroscopy study of different HTS materials with 
different orientations has been given recently by Alff et 
al. H^,^ . Details on the low temperature scanning tun- 
neling spectroscopy technique can be found in pT| , p8| and 
in the references cited therein. Typical conductance ver- 
sus voltage curves obtained by low temperature scanning 
tunneling spectroscopy are shown in Fig. M. These curves 
can directly be compared to those of Fig. |^ obtained from 
experiments using GBJs. We note that the conductance of 
the two junction types differs by several orders of magni- 
tude. Whereas one has a typical conductance of 10~^ S 
in low temperature scanning tunneling spectroscopy spec- 
tra, the conductance of GBJs ranges between 10~^ and 
10~^ S. However, this difference mainly is caused by the 
different junction area, which is by several orders of magni- 
tude smaller for low temperature scanning tunneling spec- 
troscopy. The low tenrperature scanning tunneling spec- 
troscopy data of Fig. R were obtained by performing ex- 
periments on (110) oriented surfaces of NCCO and YBCO. 
With respect to the zero bias conductance peak, the same 
result is obtained as for the GBJs. Again, only YBCO 
shows a zero bias conductance peak while such peak is 
completely absent for NCCO. We note that by low tem- 
perature scanning tunneling spectroscopy also (100) and 
(DDI) oriented surfaces have been probed. Also for these 
surfaces a zero bias conductance peak never has been ob- 
served for NCCO However, for YBCO in many cases 
a zero bias conductance peak could be detected for (100) 
and (001) oriented surfaces. For a (100) oriented surface, 
where in the ideal case no zero bias conductance peak is 
expected, this is most likely related to th e finite roughness 
of the surface as discussed in section 2.3. This effect is sim- 



ilar to the observation of zero bias conductance peaks for 
a-axis or even c-axis oriented planar junctions. 

We also would like to address another detail of the 
conductance versus voltage curves that are similar for 
both the SIN junctions used in low temperature scan- 
ning tunneling spectroscopy and the SIS-type GBJs. For 
NCCO the G{V) curves show a more pronounced super- 
conducting gap structure with a density of states that is 
reduced by about 50% within the gap. In contrast, for 
YBCO this reduction is smaller and (if at all) at the gap 
voltage only very faint and broad conductance peaks are 
observed. Finally, we note that in contrast to the spec- 
tra of the GBJs the low temperature scanning tunneling 
spectroscopy spectra show some asymmetry that is most 
pronounced in the case of YBCO. This feature has been 
observed in many spectra usinglow temperature scan- 
ning tunneling spectroscopy p5| , p9[ . So far, it is not clear 
whether this effect is related to the order parameter sym- 
metry. However, one is led to think in this direction be- 
cause the asymmetry is absent for NCCO. 
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Fig. 7. Low temperature scanning tunneling spectroscopy 
spectra for (110) oriented surfaces of NCCO and YBCO at 
4.2 K. The data are taken from reference foql. 



3.5 Anomalous Meissner currents 

GBJs formed by different HTS have been successfully used 
to determine the relative change 

Aah(O) Xab{0) 

of the in-plane London penetration depth Xab with the 
high precision of only a few A||67|]. In this technique the 
shift of the side- lobes of the Ic{H) dependence of HTS- 
GBJs is measured as a function of temperature. Details 
of this precision measurement technique have been de- 
scribed elsewhere ]67[ |. Recently, an anomalous temper- 
ature dependence of AXab was measured for a symmetric 
36° [001] tilt YBCO-GBJ |lO^. As can be seen in Fig. i, 
for temperatures below about 10 K the AXab{T) / Xab(y) 
dependence shows a clear upturn with decreasing tem- 
perature in contradiction to the behaviour expected for a 
d-wave superconductor. We note, that this anomaly can 
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Fig. 8. [Aab(r)-Aa(,(0)]/Aab(0) for a 36.8° [001] tilt symmetric 
YBCO-GBJ. The solid line shows the dependence expected for 
a superconductor with a d2,2_y2-symmetry of the order param- 
eter. The measured upturn at low temperatures most likely is 
due to anomalous Meissner currents. 
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be explained neither by assuming a pure d-wave nor a 
pure s-wave symmetry of the order parameter. In both 
cases AXab{T) / Xab{0) is expected to monotonously de- 
crease wit h de creasing temperature. However, as discussed 
in section 2.3 Andreev bound states can result in an anoma- 



lous Meissner current that is flowing in opposite of the 
usual Meissner current thereby increasing the measured 
penetration depth. This in turn results in an upturn of 
the temperature dependence of the London penetration 
depth at low temperatures. We emphasize, however, that 
a clear upturn as shown in Fig. ^ has not been observed 
for all investigated GBJs. At present it is not known what 
determines the strength of the observed upturn. However, 
we suppose that interface roughness is an important fac- 
tor. Additional experiments are required to clarify this 
issue. 



4 Discussion 

4.1 Magnetic impurity scenario 

Initially, the Anderson- Appelbaum model was used 

by most groups to give an interpretation of the observed 
zero bias conductan ce p eak in tunneling spectroscopy ex- 
periments on HTS |l||^,|l}|3|,|§l. However, later on it 
became clear that the Appelbaum- Anderson model can 
explain only a few features of the experimental observa- 
tions. For example, the predicted logarithmic dependence 
of the conductance on voltage as given by Eq. (|l|) has been 
observed for exarnple in Ref. |Q. We note, however, that 
according to Eq. {M , thermal smearing prevents the obser- 
vation of a logarithmic dependence for eV < ksT, i.e. for 
V smaller than about 1 mV at a measuring temperature of 
a few K. Therefore, the logarithmic dependence was veri- 
fied only for a relatively small voltage/energy scale. Hence, 
even though the measured data are in fair agreement with 
the Appelbaum- Anderson prediction, this cannot be con- 
sidered a definite proof for the validity of this model for the 
description of the observed zero bias conductance peak. 

We now briefly summarize the main arguments against 
the validity of the Appelbaum- Anderson model. Firstly, a 
strong argument against the applicability of the Appel- 
baum- Anderson model is the clear correlation of the tem- 
perature for which the zero bias conductance peak disap- 
pears and the critical temperature of the investigated HTS 
that varies between 20 and 90 K. Within the Appelbaum- 
Anderson model, no such correlation is expected, since 
the quasiparticle scattering off magnetic impurities lead- 
ing to the zero bias conductance peak within this model 
is independent of the onset of superconductivity. In the 
Appelbaum- Anderson model there is no distinct temper- 
ature where the zero bias conductance peak is supposed 
to disappear. Secondly, the Appelbaum- Anderson model 
cannot account for the fact that for NCCO a zero bias 
conductance peak never could be observed. If Cu^^-ions 
at the grain boundary interface are supposed to act as 
magnetic impurities causing the zero bias conductance 
peak, then it is very difficult to explain that the zero bias 
conductance peak appears only for YBCO, BSCCO, and 



LSCO, whereas it is absent for NCCO that contains the 
same copper-oxygen layers as the basic structural element 
responsible for superconductivity. We also note that the 
presence of magnetic scatterers at the grain boundary in- 
terface is not related to the presence of a supercondu ctin g 
region with a magnetic impurity concentration n,„ [101|. 
This would include a significant reduction of the critical 
temperature Tc that is clearly not observed for all of the 
GBJs. 

We next discuss the behaviour in the presence of an 
applied magnetic field. Provided that there are magnetic 
impurities at the grain boundary, an applied magnetic 
field is expected to result in a Zeeman splitting of the 
impurity levels of 26 = 2g^BH. This has been observed 
for example in the case of Ta-Ta02-Al junctions 1102]. In 



our GBJ experiments an indirect splitting of the zero bias 
conductance peak in an applied magnetic field has been 
observed. Here, indirect means that a splitting could be 
seen only in the G{V, H) — GiV, 0) curve, i.e. after sub- 
tracting the zero field curve. However, firstly at low fields 
this splitting is much too large than expected according 
to the Appelbaum- Anderson model and secondly, at larger 
fields a nonlinear increase of the splitting with increasing 
applied field is observed. In order to model this behaviour 
within the Appelbaum-Anderson model a strongly mag- 
netic field dependent (7-factor would be required that at 
low fields is an order of magnitude larger than that ob- 
served for metal - metal oxide - metal junctions. For the 
latter, typically 5 ~ 1 — 3 independent of the applied mag- 
netic field is obtained. Moreover, in Ref. jl3| a finite split- 
ting of the zero bias conductance peak was observed for a 
SIN-j unction already in zero field. Summarizing our dis- 
cussion we clearly can state that the behaviour of the zero 
bias conductance peak in the presence of an applied mag- 
netic field cannot be explained within the Appelbaum- 
Anderson model. 

In summary, for all spectroscopic experiments on HTS 
junctions reporting the observation of zero bias conduc- 
tance peaks there are strong arguments against the ap- 
plicability of the Appelbaum-Anderson model This 
statement holds for all investigated junction types with 
at least one HTS electrode, i.e. HTS/I/N-type junctions 
used in low temperature scanning tunneling spectroscopy, 
planar HTS/I/N junctions, and HTS/I/HTS type GBJs. 



4.2 BTK-model 



As discused in section 2.4, for small values of Z the BTK 



model predicts a peak in the normalized differential con- 
ductance at zero voltage that reaches a maximum value of 
2. However, this peak should decrease monotonously with 
increasing voltage and no gap structure is expected to be 
observed. That is, within the BTK-model the observation 
of both a gap like structure and a zero bias conductance 
peak is not expected. However, in all the GBJ experi- 
ments a clear gap structure and a pronounced zero bias 
conductance peak always is observed at the same time. 
Furthermore, in the GBJs no subharmonic gap structures 
at voltages 2Ajn were observed that are predicted by the 
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BTK-model. Finally, the BTK-model does not account for 
the clear differences between the s-wave material NCCO 
and the d-wave materials YBCO, BSCCO, and LSCO. 
This is evident because the BTK-model was developed 
for conventional BCS superconductors and therefore can- 
not account for effects related to d-wave superconductiv- 
ity. It has been pointed out above that the Andreev bound 
state model represents a generalized BTK-type model that 
accounts for an arbitrary Z and an anisotropic pair po- 
tential. As discussed in the next section this model well 
describes the experimental data. 



4.3 Andreev bound states in d-wave superconductors 

In the following we interpret the tunneling experiments 
on HTS within the Andreev bound state model. We will 
see that the two essential experimental findings for the 
HTS-GBJs are naturally explained by this model. Firstly, 
the clear correlation of the appearance of the zero bias 
conductance peak in the G{y) curves with the critical 
temperature is a strong hint that the presence of the zero 
bias conductance peak is related to the presence of su- 
perconductivity. Secondly, the observation that only for 
d-wave superconductors the zero bias conductance peak 
is observed gives strong evidence for the correlation of 
this effect with the symmetry of the order parameter. For 
NCCO, that is believed to be an s-wave superconductor, 
and also for conventional low temperature superconduc- 
tors a zero bias conductance peak could not be observed. 
This is fully consistent with the Andreev bound state 
model. In addition, the Andreev bound state model pro- 
vides a convincing way to explain almost all experimental 
observations of zero bias conductance peaks in different 
experimental situations within a single model [ p3| . 

We first discuss the temperature dependence of the 
zero bias conductance peak. However, we note that for this 
dependence so far no complete theoretical description is 
available at present. In Fig. ^ the normalized conductance 
at zero bias, G(0)/G„, of a LSCO-GBJ is plotted versus 
the reduced temperature t = T/Tc- For t > 0.3, the func- 
tional form is close to a 1/T dependence. Recently, such 
dependence has been predicted by Barash et al. . How- 
ever, these authors predict a 1/T dependence for small 
values of whereas we find such dependence for large t 
and, moreover, for a much wider temperature regime (up 
to t ~ 1) than predicted in Ref. |Q. With respect to this 
discrepancy we note that the faceting and the related an- 
gle averaging present in GBJs may have to be taken into 
account in order to explain the measured temperature de- 
pendence of the zero bias conductance peak. 

Next we discuss the magnetic field dependence of the 
zero bias conductance peak. In a recent paper Hu discusses 
the magnetic field effects that are expected for GBJs in- 
cluding the possibility that the magnetic field might be 
shielded or screened However, shielding effects can 
be definitely ruled out by our experimental results. By 
studying the magnetic field dependence of the critical cur- 
rent of symmetric HTS-GBJs, we always find a Fraunhofer 
diffraction pattern like dependencies. From the distance of 
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Fig. 9. Normalized conductance at zero bias of a symmetric 
36.8° [001] tilt LSCO-GBJ plotted vs. normalized temperature. 



the minima in this dependence we can conclude that the 
magnetic field penetrates the GBJ as expected. Further- 
more, the magnetic field penetrating the grain boundary 
is enha nced by a factor of order W/d due to flux focusing 
effects |l0|. Here, W is the width of the GBJ and d the 
thickness of the junction electrodes. 

We now discuss the fact that no direct splitting of the 
zero bias conductance peak in an applied magnetic field 
is observed in the GBJ experiments, whereas such split- 
ting is expected within the Andreev bound state model. 
In agreement to our observation no direct splitting of the 
zero bias conductance peak has been found also by Ekin 
et al. in experiments using planar SIN- type junctions [p2| . 
Furthermore, in experiments using low temperature scan- 
ning tunneling spectroscopy a shift or a broadening of 
the zero bias conductance peak has been observed [ p3[ . 
Only in some cases a clear split of the zero bias conduc- 
tance peak in an applied magnetic field including a split 
at zero applied field was found [ p^pl[ |. Hu proposed to 
explain these different observations with a different be- 
haviour of the magnetic field penetration However, 
since the penetration of the magnetic field and flux focus- 
ing effects are at least for HTS-GBJs well understood, this 
cannot be the reason for the absence of a direct splitting 
of the zero bias conductance peak for GBJs. However, this 
feature may be understood by considering the faceting of 
the grain boundary or, equivalently, the interface rough- 
ness present in the different junction types. Due to the 
faceting of the grain boundary barrier the GBJ is known 
to be formed by a parallel array of many small junctions 
on a micron and submicron scale with orientations that 
differ considerably from the nominal misorientation angle 
of the grain boundary (see Fig. H) . Together with impu- 
rity scattering this may lead to a masking of the direct 
splitting that is expected for junc tions with perfectly flat 
interfaces in an applied field |104|. 

Although the GBJ experiments show no direct split- 
ting of the zero bias conductance peak in an applied mag- 
netic field, there is a shift of spectral weight from zero to 
finite energies. This shift of spectral weight results in two 
peaks in the G{V, H) — G(y, 0) curves. Analyzing these 
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curves we can define a splitting 6 as half the peak-to-peak 
separation. In Fig. ^ the S values derived in this way are 
plotted versus H for a LSCO-GBJ at 100 niK and for a 
YBCO-GB J at 4.2 K. Clearly, 6 does not vary linearly with 
the applied field as predicted by the Appelbaum- Anderson 
model. For all investigated samples 5 increases linearly 
only at small fields below about 2 T. However, for larger 
fields S clearly deviates from a linear behaviour and tends 
to saturate at a constant value at high fields. This is in 
qualitative agreement with experimental results published 
very recently |[l3|,pl[ and theoretical predictions by Fogel- 
strom et al. |59|. We note that the splitting S by definition 
has to pass through the origin in Fig. That is, from 
our data we cannot conclude that there is a finite splitting 
at zero applied magnetic field as observed by Covington 
et al . , although this is suggested by the data shown in 
Fig. |1C|. The steep jump from zero to a finite splitting for 
small applied magnetic field is characteristic for all our 
samples. Furthermore, the magnetic field dependence of 
the indirect splitting derived from the G{V, H) — G(V, 0) 
curves of GBJs (see Fig. |l^) is consistent with the direct 
splitting observed by Covington et al. [|l^ for SIN-type 
junctions. The reason for this consistency is not obvious 
and still has to be clarified. 
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Fig. 10. Splitting 5 of the G{V, H) - G(V", H ^ Q) curves of 
a 36.8° [001] tilt LSCO-GBJ at T = ICQ mK and a 24° [001] 
tift YBCO-GBJ at T = 4.2 K vs. appfied magnetic fiefd. Also 
shown is the Appefbaum-Anderson-modef prediction for g = 2. 



We finally discuss a possible zero field splitting of the 
zero bias conductance peak. As discussed above, a broken 
time reversal symmetry state at the surface of a HTS can 
lead to a split of the zero bias conductance peak in zero 
magnetic field. Such direct splitting has not been observed 
in our experiments down to temperatures of 100 mK for 
LSCO and down to 4.2 K for YBCO. Our observations are 
in agreement to the recent findings of other groups |2^,^ . 
Although the detailed reason for the presence or absence 
of a zero field splitting of the ZBCO is not clear, one can 
speculate that the junction property preventing the ob- 
servation of a direct splitting of the zero bias conductance 
peak by an applied magnetic field also may result in the 
absence of any direct splitting in zero applied field. The 



most likely reason is the roughness of the surface of the 
junction electrodes, which in the case of GBJs is caused 
by the faceting of the grain boundary. Surface roughness 
together with impurity scattering may lead to a reduction 
of those quasiparticle trajectories that are responsible for 
the splitting of the zero bias conductance peak. Further- 
more, surface currents along a strongly faceted GBJ could 
tend to cancel each other. Certainly, this topic has to be 
clarified in more detail by future theoretical and exper- 
imental work. Summarizing our discussion we note that 
at present it is impossible to make a definitive statement 
with respect to the occurence of a state with broken time 
reversal symmetry in HTS-GBJs. 



4.4 NCCO-spectra 

In the tunneling spectra of NCCO-GBJs never a zero 
bias conductance peak has been observed. This observa- 
tion is consistent with an s-wave symmetry of the order 
parameter in NCCO in contrast to other HTS such as 
YBCO, BSCCO, or LSCO. Since there is additional ex- 
perimental evidence that NCCO has a dominant s-wave 
component of the order parameter, the question arises 
whether or not NCCO can be viewed as a classical BCS- 
superconductor. Furthermore, experiments performed on 
NCCO are highly important in order to clarify the origin 
of the specific features observed for the d-wave HTS. In 
this subsection we discuss the tunneling spectra measured 
for NCCO-GBJs and compare them to those obtained for 
the rf-wave HTS as well as to those expected for classical 
BCS-superconductors. 

Discussing the NCCO spectra shown in Figs. ^ and 0, 
one might think at first glance that the shape of the spec- 
tra resembles that of a typical spectrum of a BCS su- 
perconductor, if a finite lifetime of the quasiparticles is 
taken into account by introduci ng a smearing parame- 
ter r according to Dynes et al. |105|. However, looking 
more closely to the spectra it becomes evident that the 
measured G{V) curves cannot be described by a BCS 
type behaviour because of the large conductance at zero 
bias and the V-shaped gap structure in contrast to the 
U-shaped structure expected for conventional BCS super- 
conductors. The G(T^)-curves of NCCO-GBJs resemble by 
far more those of GBJs fabricated from the other (c?-wave) 
HTS materials. We note that a fit of the measured tun- 



neling spectra using the Dynes formula [105] is possible 
only in a few cases. Moreover, to obtain a reasonable fit 
r Aq has to be used what is difficult to justify. However, 
the Dynes expression cannot account for both the sharp 
peaks at the gap edge and the high spectral weight within 
the gap. Above we argued that NCCO is not a d-wave 
superconductor because a zero bias conductance peak is 
completely absent in the tunneling spectra. To account for 
the observed features an anisotropic s-wave order param- 
eter can be assumed. Then, qualitatively the high spectral 
weight at zero energy and the absence of the zero bias con- 
ductance peak can be explained. We finally note that an 
extended s-wave symmetry of the order parameter cannot 
be fully ruled out at present, since the detailed influence 
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of the surface geometry on the tunnehng spectra is not 
known. 



5 Conclusions 

Tunneling spectroscopy on differently oriented surfaces of 
HTS can be used to study the symmetry of the order pa- 
rameter, since the tunneling spectra are sensitive to the 
phase of the pair potential. Our comprehensive study of 
various HTS-GBJs gives clear evidence for the presence of 
Andreev bound states at the surface of the jimction elec- 
trodes. For the materials YBCO, BSCCO, and LSCO al- 
ways a zero bias conductance peak due to Andreev bound 
states was observed what is consistent with a d-wave sym- 
metry of the order parameter in these materials. In con- 
trast, NCCO never showed a zero bias conductance peak 
what is consistent with the absence of Andreev bound 
states and an anisotropic s-wave symmetry of the order 
parameter. Our results show a clear correlation between 
the presence of a zero bias conductance peak and super- 
conductivity in the junction electrodes. This gives strong 
evidence for an explanation of the zero bias conductance 
peak in terms of Andreev bound states and against an ex- 
planation in terms of magnetic impurity scattering. The 
presence of Andreev bound states at surfaces and inter- 
faces of d-wave HTS also explains the observed anomalous 
temperature dependence of the London penetration depth. 
Considering all tunneling data that is available at present 
from experiments using different tunneling configurations, 
the vast majority of the measured tunneling spectra can be 
well explained by taking into account the presence of An- 
dreev bound states at the surfac;e of HTS due to a rf-wave 
symmetry of the order parameter. However, some issues 
such as the presence of a broken time reversal symmetry 
state at surfaces of d-wave HTS are still unsettled. In con- 
trast to SIN-junctions, for HTS GBJs no direct splitting of 
the zero bias conductance peak in zero applied magnetic 
field was observed. 
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